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' Abstract. The set of continuous norm- preserving stochastic Schrodinger equations 

■ associated with the Lindblad master equation is introduced. This set is used to describe 
the localization properties of the state vector toward eigenstates of the environment 
operator. Particular focus is placed on determining the stochastic equation which 
exhibits the highest rate of localization for wide open systems. An equation having such 

■ a property is proposed in the case of a single non-hermitian environment operator. This 
result is relevant to numerical simulations of quantum trajectories where localization 
properties are used to reduce the number of basis states needed to represent the system 



Q\ • state, and thereby increase the speed of calculation. 

^— > 

Ctf 1 1. Introduction 

cr. 



A quantum system interacting with its environment can be described, in the Markovian 
approximation, by two complementary approaches. In the first and most commonly used 
|], 0, [| , the system is represented by a density operator and its evolution is described by 
a master equation. In the second a state vector represents the system and a stochastic 
Schrodinger equation describes the state evolution j|, ^ [J These two treatments are 
equivalent in the following sense: for all times an ensemble of state vectors generated by a 
stochastic Schrodinger equation reproduces, on average, the density operator generated 
by the master equation. The correspondence is not uniquely defined, in that for a single 
master equation there are many associated stochastic equations. 

Stochastic state vector equations, also called unravelings of the master equation, 
have been introduced in different contexts and with different interpretations. In the 
fundamental theory of quantum measurement, stochastic equations have been used 
to describe the general dynamical process of the state collapse into an eigenstate of 
the measured observable, i.e. localization j|, 0, f§ [1 0> 0' 03> 0- m quantum 



optics, stochastic Schrodinger equations have been used to describe the system state 
conditioned by measurement outcomes. In this context, an unraveling corresponds to 
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a specified measurement scheme, such as photon counting, heterodyne or homodyne 



detection |6|, TA\. More generally, in the field of open quantum systems, unravelings 
have been used as an efficient numerical method to solve the master equation 
!§, [R| ITS], |IJ, 0, [21], |22|, |23], |2§ [25], [2j]. 



The present work is motivated by a recent study of a quantum system in interaction 
with a thermal bath using the quantum jump (QJ) unraveling [27]. It is shown 
that, under some assumptions valid in the classical limit, the QJ trajectories, i.e., the 
realization of the stochastic process, approach a diffusive limit very similar to the one 
exhibited by the quantum state diffusion (QSD) trajectories. Since diffusion is expected 
to be a general feature associated with the emergence of classicality, a description of 
the whole set of continuous unravelings becomes important. This set is introduced in a 
unified way in the following section. We will show that each continuous unraveling can 
be characterized very simply by specifying its noise correlations. The set of continuous 
unravelings is then used to study how quantum state properties, such as localization, 
evolve when the unraveling changes. 

A very important characteristic of quantum trajectories which has both physical 
and numerical consequences is the localization of quantum states toward eigenstates of 
the environment operator. Working with the real noise unraveling (RN), Gisin [7j has 
shown that for an arbitrary hermitian Lindblad operator L the state vectors concentrate 
on the eigenspace of L. Percival [§] extended this result by giving a proper definition of 
the ensemble localization of an arbitrary operator and then providing analytical bounds 
for the rate of self-localization of hermitian and nonhermitian Lindblad operators for the 
quantum state diffusion unraveling. For a dissipative interaction, Garraway and Knight 
PBj , have presented numerical simulations of the localization process using the QJ 
and QSD unravelings. Starting from different quantum states, such as a superposition of 
two coherent states, a Fock state and a squeezed ground state, they have shown that such 
states are highly sensitive to dissipation. They also illustrated the localization process. 
Recently they have applied their results to describe the evolution of a Schrodinger cat 
state in a Kerr medium where localization competes with nonlinear effects |Kj (see also 

0)- 

For numerical simulation of open quantum systems, individual trajectories have 
proven advantageous over density operator computations. The main advantages stem 
from the fact that less space is needed to store and propagate in time a state vector than 
a density matrix. In addition, for trajectory methods one can exploit the localization 
property to reduce the number of basis states needed to represent the state vector, thus 
significantly reducing the time needed to calculate quantum trajectories. For quantum 
jump unravelings, when many Lindblad operators are present, it is well known that one 
can perform a unitary transformation to select one of the quantum jumps unravelings, 
in such a way as to minimize the number of basis states needed (see ref. [EJ for an 
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application of this property). The localization of the state vector for QSD has been 
exploited in the mixed representation of Steimle and al. [23] and the moving basis of 
Schack and al. pi ER EH1. 



In section |3| we use the set of continuous unravelings to describe localization 
properties for a single environment operator. Several well known results are recovered for 
a hermitian operator |4], In the case of a non-hermitian operator the minimal rate of 
localization introduced for the QSD unraveling || [TT], |TJ] is extended to the complete set 
showing that localization is a general feature shared by all the continuous unravelings, 
and a new unraveling is introduced. Some theoretical arguments supported by numerical 
simulations suggest that this new unraveling possesses the highest localization rate. 

In the present work we make use of the freedom of choice for the noise correlations 
to obtain the continuous unraveling which localizes the state vector the most. This 
transformation should not be confused with the unitary transformation discussed above. 
These two transformations are complementary and can be used together. At the end of 
section |3| we compare the localization properties of QSD and of our proposed unraveling. 

Finally, in section ^ we summarize our results and draw conclusions about the 
applicability of the set of continuous unravelings to the study of the quasi-classical limit 
of open quantum systems. 



2. Continuous unravelings 

We proceed following closely the derivation of the quantum state diffusion unraveling 
by Gisin and Percival ||. In this work, a general stochastic differential equation with 
a complex Wiener process is used as a starting point. The drift and noise terms are 
then specified by asking that the stochastic differential equation recovers on average the 
Lindblad master equation for the density operator p of the system 

p = -l[ff, p] + £ (L jP L) - \{L]L„ p}) , (1) 

where H is the system Hamiltonian and Lj (j = 1, . . . , J) the set of Lindblad operators 
which represent the influence of the environment. (Since the master equation ([!]) is valid 
under a Markovian approximation, all the stochastic differential equations considered 
apply only within this approximation). The other conditions needed to specify the 
drift and noise terms are that the state remains normalized and that the stochastic 
equation shares the same invariance properties under unitary transformations as the 
master equation. This last constraint is used to prove the uniqueness of QSD among 
the set of continuous unravelings. Removing the constraint of invariance under unitary 
transformations among the Lindblad operators, we obtain the set of continuous norm- 
preserving unravelings related to the master equation. 
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2.1. Derivation of the stochastic Schrddinger equations 

We start by considering a general stochastic differential equation of the following Ito 
form which gives the variation \dip) of the state vector in a time dt 



>><ft + I>i> (E^wJ (2) 

i=l \n=l / 



where |t>) and \uj) are vectors, <x, n are complex numbers and dWj n independent real 
Wiener processes |33[ which obey the following relationships 

M(dW jn ) = dW jn dW km = 5 jk 5 mn dt dW jn dt = (3) 

where M represents the ensemble average. The two conditions to be respected by the 
previous equation (|2]) are (i) the state is normalized for all times (ip\ip) t = 1 and (ii) for 
each time, the mean of the projector associated to the state \ip) gives the density matrix 
p = M(\ijj)(ifj\) with the density matrix p evolving according to the master equation in 
Lindblad form ([I]). In the following, these two conditions will be used to relate the drift 
term \v) and the stochastic term \uj) to the state as well as giving constraints on 
the complex numbers aj n . Notice that the ctj n may also depend on the state and 
on time t. 

By following closely the QSD derivation given in reference |§, we obtain the 
expression for the drift term 

\v) = - % -H\i,) -W ( L ) L i + - 2(^V,) |V> (4) 



where (Lj)^ = (ip\Lj\ip) is the expectation value of Lj for the state \ip). The drift term 
\v) is the same as that obtained in the QSD derivation, but the stochastic vectors \uj) 
become 

\uj) = ^=t==^(3 jk (L k - (L k )j,)\i()) j = l...J (5) 

yJ2n \ a jn\ k 

which differs from the QSD derivation by the introduction of the normalization factor 
(J2n \ a jn\ 2 )~ 1 ^ 2 and the set of complex numbers (3j k . The latter are arbitrary coefficients 
of a J x J unitary matrix which arises due to the freedom of choice in the linear 
combination of vectors (L k — (Lk)^)\i/j) used to express \uj). 
Finally one gets the equation for the state vector increment 

|^) = - \H\^)dt - i J] (L)Li + (L%(Lj)i, - 2{L%L J ) \^)dt 

+ Y.{L k - (L k )^ \^)dC k (6) 
k=i 
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This equation shows that all the indeterminacy due to the coefficients ctj n and to the 
unitary transformation (J3jk) can be included in the noise terms d(j which are given by 

aQk = 2^ Pik 7= ==■ {() 

J=l yEii \ a jn\ 

It can be seen easily that they have zero mean M(dQ) = and correlations 

dQdQ = 5jkdt and dQd(k = Cj^dt (8) 

where Cjk are correlation coefficients related to the unitary transformation ((3jk) and the 
noise coefficients a in in the following way 



i 

1 2 



En « 

Cjfc = X! WikCi with q = " m |2 (9) 

i=l ^n l^i' 



and Q are complex numbers inside the unit circle |cj| < 1 for all % = 1, . . . , J (see 
Appendix). Notice that although the stochastic process is completely specified by the 
numbers fijj. and a^, that specification will be non-unique in that all sets of numbers 
that yield the same correlation coefficients Cjk through equation @ will describe the 



same stochastic process [p^l - Thus an unraveling is completely specified when the 
correlation coefficients Cjk are given. This result provides a natural classification for 
continuous stochastic Schrodinger equations associated with the completely positive 
master equation ([!]). In ref |34| a similar classification is given for all (Markovian) 



positive master equations but only for two dimensional Hilbert spaces. This suggests 
that the same classification procedure can be extended to every positive master equation 
in arbitrary dimensions. 

In the absence of a unitary transformation (j3jk = djk) and with only one Wiener 
process N — 1, the complex noises dQ are necessarily of the form dQ = e^dWji. 
For = we recover the real noise unraveling and for <f>j = 7r/2 the imaginary 
noise unraveling. The present work shows not only that any phase <f>j can be chosen 
but also that the phase can be a smooth function of the state \if)) and of the time. 
The QSD unraveling is recovered taking two Wiener processes (N = 2) with dQ = 
(dWji + idWj2)/v2 for all j. In this special case the correlations Cjk vanish. 



2.2. Unitary transformation 

Let us introduce the following unitary transformation among Lindblad operators 

Lj = ^2u jk L k - \jl (10) 

k 
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where Ujk and Xj are complex numbers and (ujk) a unitary matrix |8|, |9|, |35|. With this 



transformation the noise terms become d(k — u jkdQ with the correlations 

j 

d(jdQ = S jk dt and dQdC k = ^ u mj u nk c mn dt (11) 

m,n=l 

These correlations will depend on the unitary transformation (ujk) unless all the 
correlation factors vanish, i.e. Cjk = 0. Since (J3jk) is itself a unitary transformation, a 
necessary condition for invariance under unitary transformation is given by 

cj = for all j = 1, . . . , J. (12) 

When only one Wiener process iV = 1 is present, the unitary invariance condition (O) 
implies <x,i = for all j. As a consequence there is no invariant unraveling with only one 
Wiener process. With two Wiener processes N = 2, the invariance condition becomes 
cP- x + (Xj2 — 0. The norm of the two complex numbers is the same |<x,i| = ja^l and the 
phases are related by <f>jx — <pj2 = 7r/2 + nir where n is any integer number. The simplest 
case n = leads to 

dlj = {p^p^ (13) 

which correspond to the complex noise used in the QSD unraveling when the phases 
4>j are set to zero. The simplest case which can satisfy the invariance condition (O) 
is given by the QSD unraveling. The phases (f)j and other choices of n introduce only 
irrelevant phase factors which can be neglected. This is the uniqueness result of Gisin 
and Percival for QSD. If one considers more than two Wiener processes iV > 3, it 
is possible to construct other unravelings invariant under unitary transformation. For 
instance: 

f> dW x + e i7T / 3 dW 2 + e~ i7T / 3 dW 3 , , 

dt = -r 3 (14) 

where we have omitted the index j and the phase factor. Since all these unravelings 
have the same correlations d( 2 = and \d(\ 2 = dt, they are equivalent ]33| and can be 
replaced by the QSD unraveling. 



3. Localization 



As an application of the set of continuous unraveling obtained in the present work, 
one can compute the rate of self-localization of a single environment operator L for a 
wide open system, i.e. H = 0, and determine the effect of the noise correlation on 
the rate of self-localization. The rate of self-localization is defined as the rate at which 
the ensemble average of the quantum mean square deviation decays 0. It is also the 
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ensemble average rate at which the state vector tends towards one of the (right- 
) eigenstates of the Lindblad operator. The quantum mean square deviationf of the 
operator L is defined as o~ 2 (L) = (L^L)^ — (L^)^(L)^. More generally the quantum 
covariance of two operators for the state \ip) is a(A, B) = (A^B)^ — (A ! )^(B)^ ||. Note 
that the quantum covariance of L with itself is just the quantum mean square deviation 
<r 2 (L) = a(L,L). We restrict our attention to a wide open system because we want to 
describe the localization process, independently of the action of the Hamiltonian. This, 
clearly, is only a first step towards a proper understanding of localization which should 
involve Hamiltonian effects as well. 

3.1. Hermitian environment operator 

For a wide open system with a hermitian environment operator L = the state vector 
\ip) evolves according to 

|#> = -~ {pL + {L%(L)^ - 2(L%L) \1>)dt + (L— (L)$) IVR (15) 

where d( is a noise of the kind previously described by (0) and (^) with an associated 
correlation factor c given by d( 2 = cdt. When the state evolves according to equation 
(|15|) we can compute the change in the expectation value of L 

d{L)^, = 2o- 2 (L)Re(d() (16) 

and the change in the quantum mean square deviation 

da 2 (L) = -2a 2 (L) 2 (1 + Re(c)) dt + 2a(L A L A , L)Re(d() (17) 

where we have used the notation La = L — (L)^. The diffusion of the expectation value 
{L)$ and the quantum mean square deviation is produced only by the real part of the 
noise term d(. 

Taking the ensemble mean shows that the expectation value M(L)^, = Tr(pL) 
remains constant and the quantum mean square deviation evolves as 

= -2 (1 + Re(c)) M (a 2 (L) 2 ) (18) 

The noise correlation c is a characteristic signature of the chosen unraveling. For 
c = 0, the quantum state diffusion result, giving a minimal localization rate of 2, is 
recovered M. In this case, as in almost all cases, the mean square deviation tends 
to zero, thus the state evolves towards one eigenstate of L. The real noise (RN) 
unraveling c = 1 is clearly the one which gives the highest rate of self-localization. As 
a consequence, for numerical simulations involving an arbitrary Hamiltonian and one 
hermitian environment operator, the RN unraveling should be used since it will produce 

f Note that the quantum mean square deviation is not an ensemble average. 
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the fastest localization (for continuous unravelings) . In the opposite case to the RN 
unraveling, if one uses the imaginary noise unraveling c = — 1, the mean localization 
rate vanishes and the state does not evolve towards an eigenstate of L. Recovering these 
well known results 0, |32[ confirms the validity of equation (|II)D. 



3.2. Non hermitian environment operator 

We consider the case of a single non hermitian Lindblad operator. Since this case is 
more difficult to treat, we restrict ourselves to the more specific case of an annihilation 
operator L = s/kcl. As for a hermitian operator, we want to determine which is the 
unraveling with the highest localization rate and find out if there is any unraveling which 
does not localize. 

The state vector evolves according to equation ([15]) and the change in the quantum 
mean square deviation a 2 (a) = (a^a)^ — (a')^(a)^ is given by 



da 2 (a) = - k (cr 2 (a) + a 2 (a) 2 + |cr(a f ,a)| 2 + 2a 2 (a)Re{ca(a\ a)}) dt 

+ 2y / ^Re | (a(a)a, a) — (a))^a(a\ a) — (a}^a 2 (a)^j d(} (19) 

which involves the quantum covariance o~(a\ a) = (a 2 )^ — (a) 2 . The equation for o~(a\ a) 
can also be derived to give 

do~(a\ a) = — k ((1 + 2o~ 2 (a))a(a\ a) + ca(a\ a) 2 + c*o~ 2 (a) 2 ^J dt 

+ 2v^Re{((T(a 2 ,a) - 2(a t )^<r 2 (a)) d(] (20) 

For the QSD unraveling c = 0, the two equations (19) and (20) are known to describe 
the localization of the state lib) towards a coherent state in the case of a harmonic 



oscillator. Furthermore, this localization is known to be globally stable [|12 |. 

Taking the ensemble mean over equation (19) removes the noise terms but introduces 
statistical correlations since the statistical mean of a product is in general different from 
the product of the statistical means. Thus one cannot obtain an immediate result for 
the mean rate of localization. However, one can notice that the drift part of da 2 (a) can 
be written as a sum of positive terms: 

a 2 (a) + (a 2 (a) - \a(ala)\) 2 + 2a 2 (a)\a(a\ a)\ ( 1 + Re{c p^ % }) (21) 

V \a(a\a)\ J 

the third term of this expression being positive since |c| < 1. As a consequence, the 
argument for global stability of coherent states, 

da 2 (a) o, . , 

M 1 ; < -KMa 2 (a), (22) 

(JjL 

is valid for all the continuous unravelings. A localization rate of k, associated with 
the exponential decay of quantum correlations Ma 2 (a) ~ exp(— Kt), is expected to be 
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independent of the choice of unraveling. This result shows that, for an annihilation 
operator, all unravelings localize and k provides a minimal bound, independent of the 
unraveling, for the ensemble mean localization rate. 

For a hermitian operator, the unraveling which localizes the most was easy to 
find since the evolution of the quantum mean square deviation is not coupled to any 
other moment. Furthermore the correlation factor c factorises, making the unraveling 
independent of the state. In the present case the situation is more complex, since 
none of these two simplifying conditions are satisfied. In the case of an annihilation 
operator, we adopt the following technique. Instead of considering the localization of an 
arbitrary state we restrict our attention to squeezed states. We will show that every 
unraveling fll5| ) with L = y/Ra preserves the set of squeezed states, i.e., if the initial 
state is a squeezed state it will evolve into a squeezed state. This has been shown for the 
QSD unraveling in ref. |12 |. The unraveling which reduces the squeezing most efficiently 
is determined. Finally, some arguments will be given as to why this unraveling should 
be the one with the highest localization rate for an arbitrary initial state. 

Squeezed states are defined as the states |7 t , a t ) which satisfy the relation 

(a - 7 t a f - a t )\j t ,a t ) = (23) 



where 7 t and a t are complex numbers which label the squeezed state |7 t , a t ) f36|. When 
the squeezing parameter j t vanishes, squeezed states reduce to coherent states. Amongst 
the properties of squeezed states, we will use for our present purposes only the relations 
between the mean square deviation and the squeezing parameter 
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v 2 s (a)= - m , 9 = lt a s (a\ay (24) 



where the index s specifies that the mean square deviation is taken with respect to 
the squeezed state \^ t) Oi t ). This last relation tells us that the mean square deviation 
depends only on the squeezing parameter. 

A condition to check if squeezed states are preserved can be obtained by 
differentiating (|23"D [jpj] , p5[| . In order to simplify the calculation, the Stratonovich 
formalism is used. In this formalism, the usual differentiation rules apply. Thus from 
(p3|), a state \ip) initially squeezed will remain squeezed if it is possible to write 

(a — 7 t a^ — a t )\dip) = {d^ t <^ + da t )\il)) (25) 



where \dip) is to be expressed in Stratonovich form. From equation ( fTSj ) and using the 
usual conversion formula X o dY = XdY + ^dXdY between the Stratonovich and Ito 
formalism [37], one can express the differential increment for the state vector as 

|#> = - \ {L ] L - 2<L%(L - (L) f )} ty)dt 

-UL"-(L%-2(L)4L-(L) lp )}\^dt 
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+ — (L)$)\ij)) o d( (26) 

Inserting this expression in the condition one finds not only that squeezed states 
are preserved but also the equations for the squeezing parameters are 

d<y t = - «7 t (l + cry t )dt (27) 

da t = - -a t dt + «7t(a t ) s (l + C7t)dt + v^Tt^C (28) 

written in Ito form. Since the evolution of the squeezing parameter ^ is deterministic, 
it is easy to find the unraveling which produces the fastest squeezing decay. It is given 
by the following correlation factor 

c = ^- (29) 
\Jt\ 

Notice that this unraveling depends on the state itself. In order to produce the maximal 
decay in squeezing, the noise term in the stochastic Schrodinger equation has to evolve 
according to the prescription given in (|29|). For such an unraveling, the norm of the 
correlation factor c is kept constant at its maximal value |c| = 1 and the phase varied 
in time along each single trajectory in order to minimize the squeezing. 

For all the unravelings with a correlation factor having such phase but an arbitrary 
norm {c r = 7"y * / 1 Tt I ? r £ [0, 1]}, the localization rate can be easily computed. This set 
includes QSD which as r = and the adaptive unraveling (p9|) which has r — 1. For 
every unraveling in this set, the squeezing parameter j t decays as j t ~ e~ Kt for large 
enough times. Using the relation ([24]), one obtains 

a 2 (a)~ e - 2Kt a(a\a)~e- Kt for t > re* (30) 

Compared to the minimal localization rate of k, this result shows that the unravelings 
which have a correlation factor c with the proper phase, i.e. all unravelings in the above 
mentioned set, produce a localization rate on squeezed states which is twice the minimal 
localization rate. Furthermore, the quantum mean square deviation a(a\a) decays at 
the same rate as the energy. Thus, the time needed for a squeezed state to become a 
coherent state is of the same order as the dissipation time, making the squeezed state a 
favourite basis for numerical simulations. 

Squeezed states are frequently represented in phase space using distribution 
functions. The Q distribution of a squeezed state is a Gaussian whose contour forms an 
ellipse. It can be easily seen that the phase fl29|) of the correlation factor c is equal to 
e «20 w h ere j s th e angle between the real axis and the major axis of the ellipse. The 
correlation factor c is defined, irrespective of its norm, as the square of the noise term 
d( (see (§)). Thus the phase of the noise term d(, which produces the fastest squeezing 
decay, is the same as the phase of the major axis of the ellipse, this later being defined 
as e 1 ^. One can say that the unraveling varies in time in order to distribute the noise 
fluctuations along the axis of largest spread of the squeezed state. 
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If the system state is not a squeezed state, the previous derivation does not apply 
anymore. Nevertheless, we can try to generalize the result for an arbitrary state. Using 
the relation (2J) between squeezing parameter and mean quantum deviation, the same 



unraveling can be specified as 

(31) 



er(a' , a) 



|cx(at, a) | 



This expression allows us to generalize the criteria obtained for a squeezed state to an 
arbitrary state. This choice corresponds to fluctuations distributed along the major 
axis of the ellipse which will give the best approximation to the state distribution 



in phase space. Another interesting aspect of the unraveling ( pTj ) is that it is also 
the one which maximizes the drift term (|2l|) in the da 2 (a) equation. This is not a 
sufficient condition to insure this unraveling localizes the most since this condition 
neglects statistical correlations. As soon as statistical correlations can be neglected, 
this unraveling will produce the highest localization rate. 

These two arguments show that for distributions which can be well approximated 



by ellipses, the unraveling (29) will give the highest localization rate. What happens 
with other states? To answer this question, numerical simulations have been employed. 
Using three different initial states, and three different unravelings, the evolution in time 
of the quantum mean square deviation a 2 (a) has been computed. The ensemble average 
Ma 2 (a) taken over 1000 trajectories is represented in figures [L]-|3|. The initial states 
chosen are: a Fock state |24), a superposition of Fock states (|23) + |25))/\/2 and a 
superposition of coherent states (\a) + | — a))/y/2 with a = 4. The unravelings chosen 
are: the unraveling with the proper phase (|3"T|), QSD with c = and real noise with 
c= 1. 

All our simulations confirm that such quantum states are very unstable under the 
effect of dissipation |2j|, this being independent of the chosen unraveling. After a small 
time Kt < 1, the state becomes a squeezed state to a very good approximation. During 
the transition from an arbitrary initial state to an almost squeezed state, the rate of 
localization is very high. 

Once the state has become a squeezed state, the mean rate of localization is 
approximately constant for all the unravelings. The localization rate always lies in 
between the minimal rate k and the squeezed state rate 2k (see figure 4). The latter 
occurs only for the set of unravelings with the appropriate phase. For the unravelings 
with a different choice of phase, the term cjt in equation (27) will give a negative 
contribution and produce a lower localization rate. In all numerical simulations, the 
new unraveling produced the highest localization. 

Finally, one can try to generalize the previous result to an arbitrary Lindblad 
operator. In this case, the unraveling which localize the most will be the one with 
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Figure 1. Ensemble average of the quantum mean square deviation cr 2 (a) showing 
the short time scale localization. The initial state is the Fock state |24). Each curve 

represent a different unraveling: the unraveling (|3l|)( ), QSD ( ) and Real 

Noise ( — • — ). The ensemble average is computed using 1000 trajectories. The errors 
bars indicate the 95%-confidence intervals. 
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Figure 2. As figure [j], but with the initial state in a superposition of two Fock states 
2- 1 /2(|23) + |25)). 



13 



20 




1 1 1 1 1 1 1 

0.01 0.02 0.03 0.04 0.05 0.06 0.07 

Kt 



Figure 3. As figure but with the initial state in a superposition of two coherent 
states 2- 1 / 2 (|a) + | - a}) with a = 4. 
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Figure 4. As figure but for a longer time scale. 
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a correlation factor given by 
\a{L\L)\ 

For a hermitian operator q, this expression reproduces the result previously obtained, 
namely 

c ~ \a( q u)\~ k 2 (<z)l { } 

for the real noise unraveling. For the annihilation operator L = a equation (|32| ) is 
equivalent to (fH|). 

3.3. Numerical simulations 

Instead of trying to justify further the unraveling ( |3"i2| ) as the most localizing one for an 
arbitrary Lindblad operator we look for an estimate of the improvement produced by 
using such an unraveling in numerical simulations. The test consists of computing 
the quantum mean square deviation cr 2 (a) along a quantum trajectory with both 
the unraveling ( |3~ll) and the QSD unraveling. The ensemble average Ma 2 (a) which 
gives a measure of the size of the wave packet as well as the variance Var{cr 2 (a)} = 
M(a 2 (a) - Ma 2 (a)) which measures the fluctuations of the wave packet size are then 
compared between the two unravelings. The quantity a 2 (a) is taken here as a measure 
of the size of the wave packet since it corresponds formally to the excitation number 
of the state displaced at the origin ||27|| . For a number states basis, this measure is 
proportional to the number of basis state needed to represent the state The system 
we used for this test is a kicked anharmonic oscillator 

H = ih(3(t)(a) - a) + ~%x<J 2 a 2 (34) 

subject to dissipation L = ^fka. The driving (3{t) is a periodic sequence of rectangular 
pulses of height /3o, length t\ and separation T2- In this case the Hamiltonian can also 
play an important role in the localization process. This system has been considered in 
previous publications |38|, 39, [35| as a simple example of a chaotic system. In conjunction 



with this system, the scaling transformations: t — Xt, k — k/X, (3 = j3 , x — x/^ 3 have 
been introduced. This scaling of the parameters does not affect the classical dynamics 
but when it is introduced in the quantum equations of motion, it allows one to perform 
the quantum-classical transition, by varying A from A = 1 (quantum) to A — > oo 
(classical). Here we use this transformation for a different purpose. In the quantum 
limit, nonlinear effects produce a strong derealization which can not be compensated 
by the localization effect of the dissipation. In the classical limit, on the contrary, 
localization is expected to be dominant. Thus varying the scaling factor allows us to 
tune the relative strength of the derealization. 
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Figure 5. Ensemble average (a) and variance (b) of the quantum mean square 
deviation c 2 (a) versus the scaling parameter A for the kicked anharmonic oscillator. 
The upper curve corresponds to the QSD unraveling and the lower curve to the time 
dependent unraveling (|3l]). The errors bars take into account statistical as well as 
numerical accuracy uncertainties. 



In figure [5], the ensemble average and variance of the quantum square deviation a 2 (a) 
are represented for different values of the scaling parameter A. The values represented 
are stationary results obtained by integrating the equations of motions over typically 
2500 periods and taking the mean over a single trajectory. Such a long integration in 
time is necessary in order to obtain a proper average over the strange attractor of the 
chaotic system. The system parameters are set to the following values x — 1; A) = 2, 
k = 0.5, 7~i = 0.98 and T2 = 1 for which the classical system is known to be chaotic. The 
precision of the numerical results is estimated by repeating the calculations for different 
numerical parameters such as the time step size. 

Figure ||(a) shows for both unravelings a slowly decaying ensemble average Ma 2 (a) 
for an increasing value of the scaling A. Notice that the amplitude of motion rescales as 
A thus the ratio a 2 {a)/M(a)a) tends towards zero when A — > oo, providing a numerical 
justification for the emergence of the classical attractor. Furthermore, fig |5|(a) shows 
that the unraveling (31) reduces, compared to the QSD unraveling, the stationary value 
of the mean size of the wave packet. The reduction can be up to 20% depending on the 
scale parameter A, the largest reduction being achieved in the quantum regime. 

More important is the reduction of the size of the fluctuations shown in fig |5](b). 
The picture suggested is that each time the wave packet deviates from a coherent state 
the QSD unraveling tends to restore the shape by applying a homogeneous noise, while 
the unraveling fl3T|) adapts by applying a non-homogeneous noise in the direction of the 
largest deviation. This adaptability does not produce an important reduction of the 
wave packet size but can stabilize the wave packet more efficiently as compared to QSD. 



16 



4. Discussion 



We have introduced the set of continuous unravelings which recovers in mean the master 
equation in Lindblad form and preserves the norm of the state vector. The quantum 
state diffusion unraveling is a member of this set, being the simplest which preserves 
the same invariance properties under unitary transformations as the density matrix. We 
have seen that each single unraveling can be specified very simply by the choice of the 
noise correlations thus providing a natural classification. For theoretical purposes, it is 
useful to work with the full set of continuous unravelings since it allows one to study 
how quantities which depend on the choice of the unraveling are sensitive to this choice. 

As a first application, we have studied the localization properties when only a 
single Lindblad operator is present. In the case of a hermitian operator, the highest 
localization rate of the real noise unraveling as well as the absence of localization 
of the imaginary noise unraveling have been recovered and explained in a consistent 
way. For a non hermitian operator, namely the annihilation operator, a new time 
dependent unraveling has been introduced. It is shown analytically that this unraveling 
provides the highest localization rate for squeezed states and numerically that this 
property is also valid for more complex quantum states. This unraveling maximizes 
the localization by continuously adjusting the phase noise according to the shape of 
the wave packet. This study provides a better understanding of the localization. For 
instance, the QSD unraveling is known to have good localization properties due to its 
invariance corresponding, in some sense, to a homogeneous distribution of noise. We 
have seen that the localization rate can be increased by maximizing the norm of the 
noise correlation factor and adjusting continuously its phase, this last constraint leading 
to the introduction of a time dependent unraveling. 

Since the new unraveling increases localization it is a good candidate for numerical 
simulations of quantum trajectories and for the solution of the master equation. A 
numerical comparison of the wave packet size and fluctuations between QSD and the 
new unraveling shows that the new unraveling performs better than QSD by stabilizing 
the size of the wave packet. 

In connection with the study of the quantum-classical transition, a recent work by 
Brun et al [27| has shown that the Quantum Jump unraveling tends to a continuous 
unraveling. It can be easily seen that this unraveling is a member of the set introduced 
in the present paper. We have shown that for a simple quantum system subject to 
dissipation all members of the set of continuous unravelings localize with a minimal 
rate given by the dissipation rate, making localization a general property valid for all 
unravelings instead of only some particular ones. 
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Appendix A. Properties of the noise correlations 

In the case of a linear combination of two Wiener process N = 2, the noise term d( 
is specified by the two complex numbers a.\ and a<i which we write as ol\ = pie 1 ^ 1 and 
«2 = P2e^ 2 . The noise correlation factor becomes 



E„l"n| 2 Pl+P2 

Using R = P2/P1 and 9 = 2(0 2 — <fii), this complex number can be rewritten as 

1 + R 2 e w 



l + R 2 



If 0i = and R are kept constant and 6 is varied, the denominator will give a circle in 
the complex plane, centered at (1,0) and of radius R 2 . The denominator will restrict 
the circle to be inside the unit disk, centered at the origin. The phase 0i produces only 
a rotation around the origin. Thus the number c can take any value inside the unit 
circle only: \c\ < 1. Furthermore, this result can be generalized without difficulty for 
an arbitrary number N of noise terms. 
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